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Opto- and electromechanical systems offer an effective platform to test quantum theory and its 
predictions at macroscopic scales. To date, all experiments presuppose the validity of quantum me¬ 
chanics, but could in principle be described by a hypothetical local statistical theory. Here we suggest 
a Bell test using the electromechanical Einstein-Podolski-Rosen entangled state recently generated 
by Palomaki et al. [1] , which would rule out any local and realistic explanation of the measured data 
without assuming the validity of quantum mechanics at macroscopic scales. It additionally provides 
a device-independent way to verify electromechanical entanglement. The parameter regime required 
for our scheme has been demonstrated or is within reach of current experiments. 


Introduction .—The interaction of light with a mechan¬ 
ically compliant mirror was at the heart of a number of 
Gedankenexperiments in the early days of quantum the¬ 
ory, and still represents a textbook paradigm illustrating 
the basic principles and intrinsic limitations of the mea¬ 
surement process in quantum mechanics [2, 3]. In recent 
years experiments in the field of opto- and electrome¬ 
chanics managed to approach some of these textbook 
examples of controlled quantum dynamics to an amaz¬ 
ing degree: Sideband cooling of mechanical oscillators 
close to the ground state [4, 5], measurement back-action 
noise [6, 7], ponderomotive squeezing of light [8, 9], co¬ 
herent quantum state transfer [10, 11], feedback control 
within the thermal decoherence time [12], and the genera¬ 
tion of Einstein-Podolski-Rosen entangled states [1] have 
all been realized with nano- to micron-sized mechanical 
oscillators coupled to optical or microwave fields. All 
of these experiments impressively demonstrate the prin¬ 
ciples of quantum mechanics working at astonishingly 
macroscopic scales. 

Although all of these effects perfectly confirm the pre¬ 
dictions of quantum theory, it is important to note that 
the data taken in the corresponding experiments could 
be perfectly well explained in terms of a hypothetical 
local statistical theory, that is, by a theory assuming lo¬ 
cal hidden variables that fundamentally determine the 
ostensibly random measurement results observed in an 
actual experiment. The possibility of such an alterna¬ 
tive, deterministic and local explanation can be ruled 
out by performing specially designed tests such as the 
violation of a Bell inequality [13]. In their simplest form 
Bell inequalities present constraints on correlation func¬ 
tions from measurements on the two spatially-separated 
halves of a bipartite system, which result from the very 
assumption of local hidden variables. By making mea¬ 
surements on an entangled quantum system a violation 
can be achieved, thereby excluding any explanation of the 
observed data based on a local hidden variable theory. 
Vice versa, the violation of a Bell inequality guarantees 


the existence of entanglement in a measurement-device- 
independent manner [14]. 

Here we propose a scheme to violate a Bell inequality 
using an electromechanical system, which enables us to 
test local hidden variable models at macroscopic scales. 
The violation can be achieved with the Einstein-Podolski- 
Rosen entangled state recently reported in [1] among 
a mechanically compliant capacitor and a microwave 
pulse. In this electromechanical setup high-efHciency 
non-number resolving detectors for photons and, indi¬ 
rectly, also for phonons can be realized by coupling to a 
superconducting qubit, as illustrated below. These tools 
suffice to violate a Bell inequality of the Glauser-Horne- 
Shimony-Holt (CHSH) type originally introduced in [15], 
which was previously violated with two-mode squeezed 
optical fields [16]. We show that a significant Bell vio¬ 
lation can be attained with parameters that are close to 
the values from [1] taking into account the most domi¬ 
nant channels of loss and decoherence (e. g., cavity losses, 
photon loss in transmission lines, thermal decoherence of 
the mechanical oscillator) non-perturbatively. 

In view of the recent experiments confirming predic¬ 
tions of quantum theory for nano- and micromechanical 
systems we suggest, much in line with [17], to strive now 
for the next level of tests of quantum mechanics chal¬ 
lenging classical assumptions of realism and locality at 
macroscopic scales. 

In the following we will first introduce the specific type 
of Bell inequality relevant to our scheme, then we explain 
the experimental setup and the protocol for violating this 
Bell inequality, and finally we present a detailed quanti¬ 
tative model from which we infer prediction for the de¬ 
gree of Bell violation that can be expected under realistic 
conditions. 

Bell Inequality .—Let 0 ^ 1 ( 0 :) and iub{P) be a set of ob¬ 
servables for two quantum systems A and B, labeled by 
measurement settings a and /3. Each observable has 
possible measurement outcomes ±1. Under the pre¬ 
misses of realism and locality the correlations Eap = 


2 


((U^(a) 0 (Ub(/ 3)) between pairs of measurements for two 
settings ai( 2 ) and /3i(2) obey the Bell (CHSH) inequal¬ 
ity [13] 


2 ^ |-£'cti/3i + -E'ai/32 “f -E'q2/3i -^02/32 I —’ (1) 

With a suitable choice of observables <i7yi(a) and vsiP) 
measured on appropriate entangled quantum states this 
bound may be violated. The maximal violation allowed 
by quantum mechanics is S' = 2\/2 [18]. 

Here we consider measurements on continuous-variable 
systems performed with a detector that allows to distin¬ 
guish the vacuum state from all other Fock states. For 
a mode of an optical field this corresponds, for example, 
to a standard single-photon counter. If a coherent am¬ 
plitude —a is added to the mode before it is detected, 
the measurement effectively distinguishes between ele¬ 
ments of the positive-operator-valued-measure (POVM) 
{Pa,t — Pa}, where Pa = |a)(al denotes the projection 
operator onto the coherent state ja). The technique of 
measuring this POVM is commonly referred to as weak 
field homodyning [19, 20]. Let the detection of ja) cor¬ 
respond to the measurement result -1-1 and the comple¬ 
mentary event to —1; the observable is then effectively 
described by 

©■(a) = la)(aj — [l — la)(aj] = 2ja)(aj — 1. (2) 

The correlation functions can thus be expressed as Eap = 
4,{Pa®Pp)—2{Pa®t-\-\®Pij) + l. Observables of the form 
(2) and the corresponding Bell inequality (1) have been 
introduced first in debates regarding the nonlocal proper¬ 
ties of spatial superpositions of single photons [15, 21, 22], 
and have been realized experimentally in [23]. Remark¬ 
ably, the Bell inequality (1) can be violated not only with 
non-Gaussian states (such as single photon states), but 
even with Gaussian entangled states. For example, a two 
mode squeezed state 

= sechr ^ (—tanhr)" \n)A ® \n)B (3) 

n 

yields a maximal violation of 5” « 2.45 for 6.3 dB of 
squeezing (r « 0.76) for optimized values ai( 2 ) and 
i5i(2) [15j 24, 25]. An experimental demonstration with 
squeezed light was reported in [16]. 

An electromechanical two-mode squeezed state has re¬ 
cently been realized in [1]. In the following we will 
show that this EPR-entangled state shared between a 
micron-sized mechanical object and a travelling-wave mi¬ 
crowave pulse can be used to violate the Bell inequality 
(1). To perform a measurement of (2) on the electro¬ 
magnetic mode, we employ a qubit integrated into a mi¬ 
crowave cavity, which can directly be used as a single- 
photon detector. In order to achieve single-p/ionon de¬ 
tection on the other hand, the mechanical state has to 
be first transferred to the microwave field. Assuming the 
mechanical state is faithfully transferred, another qubit 
then effectively acts as a single-phonon detector. In this 
way photon-phonon correlations between a microwave 


pulse and the mechanical oscillator can be inferred from 
photon-photon correlations between two pulses. 

The Protocol .—The proposed protocol can be summa¬ 
rized as follows [see Fig. 1(a)]. We hrst generate elec¬ 
tromagnetic EPR entanglement between the mechanical 
mode and a microwave pulse (A), by driving the elec¬ 
tromechanical system on the blue sideband. The me¬ 
chanical state is then swapped to a second pulse {B) by 
employing a red-detuned drive [26]. Using two microwave 
cavities containing qubits we can subsequently measure 
the observables €A{oi), < 1 ^ 3 {(3) on the two pulses and cor¬ 
relate the measurement results. The protocol thus ef¬ 
fectively consists of three steps, which we first discuss 
in an idealized scenario. Perturbative dynamics will be 
discussed afterwards. In the following we denote by Ci, 
cj the bosonic annihilation and creation operators obey¬ 
ing [ci,c|] = 5ij (where i,j label subsystems as detailed 
below). 

(i) Entanglement Generation: The mechanical oscil¬ 
lator is pre-cooled to its quantum ground state by pas¬ 
sive sideband cooling as demonstrated in [4]. The elec¬ 
tromechanical system is then driven by a blue-sideband 
pulse for a time ti , generating the entangling, two-mode 
squeezing dynamics iJgq = <7sq(cicCm + clc'^n) between 
the mechanical mode (m) and the LG mode (Ic) [1]. For 
an electromechanical coupling psq much smaller than the 
(energy) decay rate kjc of the LG circuit (weak-coupling 
regime), entangled photons leave the cavity faster than 
they are created, such that the mechanical oscillator be¬ 
comes entangled with a travelling microwave pulse with 
an exponentially growing temporal profile oc exp(rsqt), 
where Fgq = Ag'^^/nic [27]. The mechanical oscillator and 
the light pulse will then approximately take on a two¬ 
mode squeezed state of the form (3), with r = FsqTi. 

(ii) Photodetection: An observable of type (2) can be 
measured on the microwave pulse as follows. We assume 
a superconducting qubit is integrated into a cascaded mi¬ 
crowave cavity (c). It is initialized in its ground state \g), 
and exhibits a switchable dispersive interaction [28-30] 
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FIG. 1. (a) Electromechanical circuit (formed by the LC res¬ 
onator and the mechanical oscillator) and cascaded microwave 
cavities with integrated qubits, used as single-photon detec¬ 
tors. (b) Level scheme of superconducting qubit coupled to 
the intracavity microwave field. The solid green line repre¬ 
sents the resonant transition which we seek to drive, while 
the dashed line indicates an off-resonant transition. 
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Hq = xit)(J 2 c\cc with the cavity, where is the Pauli z- 
matrix. After step (i), assuming that the cascaded cavity 
possesses a bandwidth sufficiently large to accommodate 
the exponential pulse, the pulse enters the cavity. If at 
this point the qubit-cavity interaction is switched on, the 
qubit transition frequency is shifted by xric for ric cavity 
photons. Given that the shift per photon x is larger than 
the linewidth of both the qubit and the cavity (this shift 
can approach values of 1000 times the linewidths of qubit 
and cavity [31]), then the qubit can be flipped selectively 
from its ground to its excited state by applying an exter¬ 
nal TT-pulse that is resonant within the ric = 0 subspace 
[see Fig. 1(b)]. A subsequent measurement of the qubit 
thus allows us to distinguish the vacuum from all other 
Fock states. Adding a coherent displacement of —a to 
the transmission line before the cascaded cavity (or an 
appropriate amplitude to the cavity directly), results in 
effectively measuring the observable (2). 

(hi) Phonon-Detection: In order to perform a mea¬ 
surement of (2) on the mechanical oscillator, its quan¬ 
tum state is swapped to a second microwave pulse and 
the measurement of step (ii) is repeated. The state swap 
can be achieved by driving the electromechanical sys¬ 
tem with a red-sideband pulse, which generates a beam¬ 
splitter like interaction Fibs = Sbs(c]),Cm + ClcCL) [11]- 
These dynamics create a microwave pulse whose quan¬ 
tum state ideally is identically to the state of the me¬ 
chanical oscillator at the end of step (ii) [27]. How¬ 
ever, due to its exponentially decaying temporal envelope 
oc exp(—Fbst) (with Fbs = 4(7bs/^ic)i the pulse will be ab¬ 
sorbed by the cavity rather poorly. In order to avoid the 
associated photon loss we require that both the strength 
of beam-splitter coupling ghs{t) and the cascaded cav¬ 
ity’s linewidth Kdt) can be tuned as a function of time. 
The coupling strength can be tuned by tailoring the in¬ 
tensity of the drive field incident on the LC circuit; a 
time-dependent coupling between cavities and transmis¬ 
sion lines (i. e., a time-dependent linewidth) has recently 
been demonstrated in [32-37]. For an optimized control 
sequence (see App. B for details) the resulting dynamics 
approximates an ideal quantum state transfer from the 
mechanical oscillator to the cavity. Adding an appropri¬ 
ate coherent amplitude —/3 to the microwave pulse thus 
provides an effective measurement of (2) on the phonon 
mode via a measurement performed on the microwave 
pulse. 

From the statistics obtained by repeating steps (i) to 
(iii) for fixed amplitudes a and j3 one can compute the 
correlation between the two pulses, which repre¬ 
sent photon-phonon correlations between the first mi¬ 
crowave pulse and the mechanical oscillator. Performing 
the procedure for appropriate amplitudes (measurement 
settings) ai( 2 ) and /3i(2) ultimately allows to violate the 
Bell inequality (1). 

The Model .—To show that a violation of a Bell inequal¬ 
ity can be achieved in state-of-the-art electromechanical 
experiments [1, 11] we provide a detailed model of all 
steps, including the dominant decoherence channels. In 


particular we include non-perturbatively mechanical de¬ 
coherence, photon losses, and counter-rotating terms of 
the radiation-pressure interaction. In order to model the 
measurement of (2) using one of the microwave cavities 
containing a qubit, we treat the cavity as a cascaded 
system [38, 39], to which the LC circuit couples unidirec- 
tionally. This allows us to correctly describe the transfer 
of the pulse into the cavity without treating it explicitly. 
The state of the three modes (mechanics, LC circuit, one 
of the microwave cavities) is described by the density ma¬ 
trix /i, whose evolution during steps (i) and (iii) follows 
the master equation 

fi = -i[iOmc\nCm “ AcJ^Clc -f (^Qc -b g* c\^){Cm + c\^), M] 
-b Cmh' + KlcT’[cic]ti -b KcV[Cc\pL 

- \/AtKlcKc/4 I [cj, Clc^] -b [pLc\^, Cc] | . (4) 

In the Hamiltonian dynamics (first line) LOm denotes the 
mechanical frequency, A = Wdrive — wic the detuning 
between the frequencies of the LC mode wic and the 
drive field Wdrive, and g{t) = \J D(t)Kie/2Fa;ic 

is the linearized optomechanical coupling. The cou¬ 
pling per single photon is denoted by g^ and P(t) is 
the power of the drive field which may vary slowly in 
time as long as P/P <C max(Kic, |A|) [27]. The sec¬ 
ond line describes decoherence processes by means of 
Lindblad operators 'D[a]pL = apLo^ — ^a/ag, — go/a and 
Cmg = lm{n + l)T>[cm\g + 7 mhT>[cJ,^]/r. The full-width- 
at-half-maximum damping rate of the mechanical oscil¬ 
lator is 7 m and n = [exp{fujJm/kBT) — 1]“^ is its mean 
occupation number in thermal equilibrium at tempera¬ 
ture T. The third line models the cascaded coupling 
of the electromagnetic system into the microwave cavity 
[38, 39]. The efficiency of the transmission channel is A*. 

The master equation (4) describes Gaussian dynam¬ 
ics, and can in principle be integrated exactly. In order 
to speed up integration and numerical optimization, we 
adiabatically eliminate the LC circuit (valid in the weak- 
coupling regime g k\c) and integrate the dynamics in a 
frame rotating at the mechanical frequency ujm, which is 
by far the fastest time scale in the problem (see App. C). 
In step (i) the detuning is chosen on the first blue side¬ 
band, A = Wm, yielding the effective master equation 

P — ^mP “b crsq27[Cm]p “b T)\y/P 

T (1 '^t)rsq^[Cm]p -b — y/AtFgqKc [CmCc -b C^C^, p\ 

(5) 

for p = tric(/r), with e = 1/[1 -b (dwm/Kic)^]- As it turns 
out, it is advantageous to slightly mismatch the cavity’s 
bandwidth with respect to the exponential envelop of the 
light pulse; this is due to the finite duration of the pulse, 
which causes a spectral broadening. We therefore set 
Kf. = uFgq, and optimize later with respect to v. In step 
(iii) we use a red-detuned pulse with A = —ujm, leading 
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to the equation 

P = C,mP + f^hj>[c\n]p + T>[^cCc “ i\/\^hsCm]P 

“t“ (1 ^t)^hs^\^m\p “t“ 2 -^irbs^c [^m^c “t“ P] • 

( 6 ) 

Both the linewidth Kc(t) of the cavity and the amplitude 
\/P{t) of the pulse [and therefore the effective electrome¬ 
chanical coupling strength rbs(t)] needs to be shaped as 
detailed in App. B to maximize the read-out efficiency. 

In order to evaluate the quantity S in (1) for the 
bipartite system consisting of the two light pulses, we 
in turn integrate equations (5) and (6) for durations 
Ti and T 2 respectively, assuming that initially the me¬ 
chanical mode is in a thermal state with a mean occu¬ 
pation number tiq and the respective microwave cavity 
is in the vacuum state right before the arrival of the 
pulse. As the system is Gaussian, its state is fully de¬ 
termined by the first and second moments of the vector 
X = (xm, J/mia;c,^c); where Xk and are quadrature 
operators obeying [xk,yi] = iSki- To evaluate the quan¬ 
tity S it suffices to calculate the symmetrized covariance 
matrix T,ki = \{X},Xi + XiXk) - {Xk){Xi) at the end 
of the pulse sequence. The master equations (5) and 
(6) lead to a differential Lyapunov equation of the form 
i] = FS -I- SF"'’ -|- N [40] . The explicit form of the matri¬ 
ces F and N is given in App. D. The Lyapunov equation 
is linear and can be integrated analytically [even for time- 
dependent parameters rbs(t), Kcit)]. The covariance ma¬ 
trix S determines the characteristic function from which 
the Bell inequality violation (1) can be calculated along 
the lines of [15, 25]. 

Results. —We optimize the resulting value of S with 
respect to the measurement settings ai( 2 ) and /3i(2)j the 
pulse duration ti and T 2 , and linewidth of the microwave 
cavity [in step (i) only], parameterized by v as discussed 
above. The optimization is performed for a fixed trans¬ 
mission loss 1 — At, bath temperature T and for a given 
maximal coupling ^max = g{t). To facilitate the 

comparison between different experimental platforms, it 
is instructive to parameterize this coupling strength by 
means of the cooperativity C = K\c'yTnin + ^)■ The 

results for the maximal Bell correlations are plotted in 
Fig. 2 versus C, for different values of the transmissivity 
At and the initial mechanical occupation number ng. 

We conclude that a significant violation of the Bell 
inequality (1) can be achieved with cooperativities and 
initial mechanical occupation numbers that are feasible 
in electromechanical systems. Cooperativity values of 
up to C « 300 and occupation numbers of ng = 0.34 
and 0.25 have, for example, been demonstrated in [4] 
and [41], respectively. The greatest challenge will be to 
bring the overall transfer efficiency above 90%. This is 
a lively research activity in the superconducting qubit 
community, however. 

The optimal values for ti , T 2 and v maximizing the Bell 
violation are shown in Fig. 3, from which we infer that 



FIG. 2. Bell correlations against cooperativity optimised 
with respect to ri, T 2 and v for different transmissivities 
At = 1,0.98,0.95,0.92 (red, yellow, green, blue), and initial 
mechanical occupation numbers no = 0.1,0.25 (solid, dashed 
lines). Other parameters are h = 40 and KicloJm « 1/8. 



c 

FIG. 3. Optimal values of ri (top), T 2 (middle), and v (bot¬ 
tom) for the corresponding values of S in Fig. 2. 


the optimal values of both ti and T 2 lie well below the 
effective coherence time 1 jn^m of the mechanical system. 
For increasingly long pulses (for low values of C), the 
optimal value for the decay rate of the microwave cavity 
is Kc ~ Tgq, as expected. 

Conclusion .—In this Letter we present an effective 
scheme to demonstrate the violation of Bell’s inequality 
using EPR entanglement shared between a mechanical 
oscillator and a microwave field. We analyse in detail 
the experimental implementation, including the primary 
decoherence channels, such as photon losses and thermal 
mechanical noise. We show that a significant violation 
of Bell’s inequality is achievable with electromechanical 
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systems. We want to emphasize that using an experi¬ 
mentally considerably less complex setup employing only 
a single detection setup can still be used to demonstrate 
electromechanical entanglement in a device-independent 
manner. 

Note that an equivalent scheme can be considered in 
the optical domain using conventional photodetectors in¬ 
stead of the qubit as a photon counter. During prepara¬ 
tion of this manuscript we became aware of related work 
along this line by Vivoli et al. [42]. 


Appendix A: Photodetection 


A measurement corresponding to the POVM 
{|a)(a|,l — |a)(a|} on the cavity mode c can be 
realized in the following way: Due to the dispersive 
interaction, the intracavity field shifts the qubit’s 
resonance frequency by ny, where n is the number of 
intracavity photons. We start with the qubit in the 
ground state and apply a 7r-pulse at a frequency Wgg. 
If the energy levels of the qubit are well resolved, i. e., 
the frequency shift between adjacent levels x is larger 
than the linewidth of the qubit’s excited state, we can 
describe the unitary evolution generated due to the 
TT-pulse by 


= |0,e)(0,g|,, + |0,g)(0,e|e, + ^|0(Z|c®l, (Al) 

1=1 

Measuring the state of the qubit after the 7r-pulse is then 
described by the measurement operators = {r\UTr\g)q 
for possible outcomes r G {g,e}. We thus see that mea¬ 
suring the qubit realizes a projective measurement on 
the cavity field where Me = |0)(0|c. Mg = Ic — Me cor¬ 
respond to finding the qubit in the excited and ground 
state respectively. The probability to find the outcome 
r is Pr = ti{Mi^Mrp). If before applying the 7r-pulse, we 
displace the cavity field by —a, i. e., p -G D{—a)pD^ {—a) 
and set r = e, we find 

Pe = tr{\0){0\eD^a)pD{a)) 

= tr(|a)(a|cP) = (Pa), 

which is what we need to evaluate the parameter S. 


Appendix B: Pulse Shape Optimisation 


the adiabatic Langevin equations corresponding to mas¬ 
ter equation (6). Here we are only interested in the clas¬ 
sical dynamics and we thus introduce the mean values 
/3(t) = {cm{t)) and ^{t) = {cc{t)). Their equations of 
motion are 

m = (Bia) 

iii) = -y - i\/rbs(i)«c(i)/3(i)- (Bib) 

Clearly (Bla) has the solution 

/3(t) = exp y dTrbs('r)^/3(0) (B2) 

and we can thus write, after formally integrating (Bib), 


^(^• 2 ) = e 


0-5 dsK^(s) g- 




Jo dsrbs(s) 


m- (B3) 


Our goal now is to maximize the term in brackets in the 
second line (we call it I) which quantifies the fidelity of 
the state swap and fulfills 0 < / < 1. We define two 
functions 

vit) = (B4a) 

w(0 = 

which we assume to be square integrable. We can thus 
write the overlap / as a scalar product on the underlying 
vector space, which obeys the Cauchy-Schwartz inequal¬ 
ity 

I = {v,w)'^ < {v,v){w,w). (B5) 

The right-hand side is easily evaluated and we find 

{v,v) = l-e"-^”, 

{w, w) = 1 — e~^'“ 

with Ky = dsrbs(s) and Kyj = dsKc(s). The 
inequality (B5) is saturated for the choice v = w, or 
equivalently (as v,w >0) for v'^ = w'^. A possible choice 
for Tbs and Kc is thus 

rbs(t) = Ae-l*"' 

Ke{t) = Ne- ^0 ^ (B6b) 


Here we discuss the optimal-control problem of how to 
transfer the mechanical quantum state into a cascaded 
cavity. We use a red-detuned (A = —uJm) light pulse 
of length T 2 , whose power P{t) can be varied in time 
[leading to a adiabatic coupling strength rbs(t)]. Ad¬ 
ditionally we assume we can tune the bandwidth Kc{t) 
of the cascaded cavity. To identify the relevant tempo¬ 
ral shapes we rewrite the system’s evolution in terms of 


where N fixes their norm. The set of differential equa¬ 
tions corresponding to (B6) is 

Tbs{t) = Kc{t)Tbs{t), (B7a) 

Ka{t) = -Kc(t)Tbsit), (B7b) 

with the boundary conditions 

rbs(T2) = «c(0) = N. (B8) 
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Non-singular solutions of these equations are given by 


rbs(/) = 

M 

(B9a) 

X _|_ gM(2i-T2) ’ 

Kc{t) = 

M 

(B9b) 

1 _|_ g-M(2t-T2) ’ 


where the parameter M is determined by the condition 
rbs('r 2 ) = ^c(O) = M{1 + Functions (B9) are 

shown in Fig. 4. 



FIG. 4. Optimal time-dependent valne of FbaCt) = 
45bs(t)^/Kic (yellow line) and Kc{t) (blue line) for the sec¬ 
ond pulse in terms of the mechanical FWHM linewidth 7 ™ 
(right axis) and the mechanical resonance frequency (left 
axis), where we assumed a Q-factor of Q = 3 • 10®. 


Appendix C: Adiabatic Elimination of the LC Mode 


e, which we write as JC = We can show that 

up to second order the expansion coefhcients are given 
by 

JCi{t)=VCi{t)V, (C4a) 

lC2{t)=VC2{t)V + VCi(t) j dTQCi{T)V, (C4b) 

Jo 

and /Co = 0. In order to apply this method for adiabatic 
elimination we need to identify a suitable subspace V 
which describes the relevant dynamics. This subspace 
must be chosen such that [44] 


VCo = CoV = 0 , 

(C5a) 

P£iP = 0 , 

(C5b) 

V = lim e^“*. 

t —^0 

(C5c) 


2. Elimination of the LC Mode 

To eliminate the LC mode we choose = tric(/r) ® 
Pvac which projects the state p, onto the cavity’s ground 
state. This is the subspace we are interested in, as in the 
limit g <C Atic (in this section we take g to stand for either 
5 sq or g^s) all photons scattered from the mechanical os¬ 
cillator into the resonator (typically on a timescale 1/g) 
immediately decay from it (on a much shorter timescale 
1/kic). For the microwave cavity we have Kq ~ 4:g^jK\c 
and thus Kefnic ~ (25/kc)^- It is convenient to introduce 
the parameter R = ^k\cKc- Consequently we assume a 
separation of time scales of the form 

Kc^ g,R<. Klc, Wm, Ac, Ale 


1. Time-Convolutionless Projection Operator 
Method 

Adiabatic elimination of the cavity mode on the ba¬ 
sis of the master equation (MEQ) (4) can be achieved 
through the time-convolutionless projection operator 
method [43]. Given a set of linear operators Ci {i G 
{ 0 , 1 }) and a corresponding equation of the form 

= [£0 + e£i(/) + £^£2] M(t), (Cl) 

the goal is to find the dynamics of p. on a subspace defined 
by the projection V = 1 — Q. We first transform into an 
interaction picture with £9 and define 

p(t) = exp (-£ot)/i(/), (C2a) 

£i(t) = exp (-£ot)£i(i) exp (£ 9 /). (C2b) 

The equation of motion for Vp can then be written in 
the form 

= mvm. m 

where we assumed that Qp{Q) = 0. We are only inter¬ 
ested in a expansion of /C in powers of a small parameter 


which reflects the structure of the Liouvillian if we iden¬ 
tify £ with gf Kc- We neglect the mechanical decoherence 
for now (i. e., we set 7 ^ = 0 ) and add it again in the end. 
(This can be shown to be exact.) 

By going into an interaction picture with the free 
Hamiltonian — Accjcc we can write the MEQ 

(4) in the required form p = [£9 -|-£i(/) + L 2 \p with the 
definitions 


£op = -iAic[c|'^cic, p] - 1 - Kic T>[cic]p, 

(C 6 a) 

£i(/) = e“^™* £f^p + £ff p + H. c. 

(C 6 b) 

£2 = KcT>[Cc]p, 

(C 6 c) 

and 


JJtnP = -i5[cL(cic + cfj.Ai], 

(C 6 d) 

JJfnP = (A+m)''' = -i5[Cm(cic + ej'j,/^]. 

(C 6 e) 

£ffp = -v^k[cI,cic 1 ^], 

(C 6 f) 

JJcP = {JJcPf = \/^tK[Ce, pc\^, 

(C 6 g) 

We can then show equations (C5) and additionally find 
the useful relations 

Q£oQ = 0 , 

(C7a) 

•pe^ot = g^ot-p = p.^ 

(C7b) 
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which we can use to evaluate equations (C4). We first 
introduce p{t) = tric(/r(t)) and thus have Vp,{t) = p{t) 0 
Pvac- We then immediately find /Ci(t) = VCi{t)V = 
VCi{t)V = 0 and VC 2 V = C 2 V. The second term in 
(C4b) is more involved. Taking into account that = 

0 we find the expanded expression 

VCiit) [ dTQCi{T)Vp{t) = 

Jo 

f dre'^-^rC+e^-^C-Vpit) 

Jo 

+ f dr VC-C-rp{t)e-^'--* 

Jo 

Jo 

+ [ dTe-'‘^™^iP/:+e^“^/:+^p(t)e-‘(^‘=-“’")‘+H.c. 

Jo 

The first term results from the optomechanical interac¬ 
tion alone and will give rise to the familiar heating and 
cooling terms derived in [45]. The second term is fast 
oscillating and will later be dropped in a rotating-wave- 
approximation. The last two terms correspond to scat¬ 
tering of photons into the cascaded cavity mode. Note 
that depending on the detuning of the second cavity only 
one of the mechanical sidebands is resonantly scattered 
into the second cavity mode. We evaluate the impor¬ 
tant terms separately, treating the integrand only for the 
moment. The first term gives 

trie(7^/:+= 

-g^ (tric(cic(r)c]'^Po)[cL,Cr„M]-tric(cicc|'^{T)po)[cL,MCm])- 
The third and fourth term give respectively 

tric(T’£;!‘ e^“^ £“T’/9(f)) = tric(cic(T)cJ'^po)[cJ, c^/r], 

tric(iP£+ e^«^/:+iPp(t)) = ig^/\tRtY^^{clc{T)c\^po)[cl,c\^p\■ 

As we will drop the second term in a rotating-wave 
approximation (RWA) we do not state the explicit ex¬ 
pression here. Next we need to find the correlation 
functions tric(cic(<)cj’^(0)/9vac) and tric(cic(0)c]’^(f)pvac) = 
ti'ic(cic(i)cJj.(0)yOvac)*- These can be calculated using the 
quantum regression theorem. We find 

tric(cic(t)c}'^(0)pvac) = exp [(iAic - Kic/2)t], (C8a) 

tric(cic(0)c|^(t)pvac) = exp [-(iAic -f Kic/2)t]. (C8b) 

If we neglect terms of the form exp(—Kict) we thus find 
dTe=^'“™’' tric(cic(T)c]'^(0)pvac) « (C9a) 

where r]± = — i(—Aic ± Taking this all to¬ 



gether we neglect fast rotating terms in a RWA, and find 


p = -i[SLUmclnCm, pj+T. T)[Cm]p+T + T)[cljp+Kc T’lCcj/O 
- g\/>^tKc/K\c{[cl, S+p] -HH. c.} 

- g\/>^tKc/Kic{[cl, s-p] _pH. c.}, (CIO) 


where we introduced jump operators s+ = —iKicP+cl^ 
and s_ = —iKicP-Cm- The other parameters are given 
by 6 ujm = 5 ^Im(p_ -|- p+), r± = 2g‘^Reg±. To get a 
time-independent MEQ we first go into an interaction 
picture with SuJmclnCm, which eliminates the first term 
in (CIO) and leads to the replacement —)■ e^'“™ * 

(with = U!m + Sujm)- We then introduce an averaged 
density operator, formally defined by 


p{t) dt = 




drp(r). 


(Cll) 


6 t has to be chosen to be long on the mechanical oscil¬ 
lation frequency, but short on all other timescale s. This 
leads to the inequality » 1/St » g^/nic, g^fn/jK^. 
If we now consider the cases Ac = dzuj^ and neglect 
fast rotating terms [which give a correction of order 
0{l/oj///St)] we eventually obtain for p the coarse-grained 
MEQ 


p = r_ V[Cm\p + V+V[c\^\p + Kc T>[Cc\p 

- g^/Xpi/J^{[cl,s±p] + [psl.,Cc])■ (C12) 

We can now use the identity 

I?[a-|-&]p = 'D[a]p+'D[b]p+-[a^b—ab\p] — {[a\bp] + [pb\a]) 

to convert this to explicit Lindblad form. For the two 
cases Ac = this leads to 

p = CmP + i\/AtKcr+/4[e““^+ CmCc + H. c., p] -f r_ 'D[cm\p 

+ (1 - At)r+ X>[cj„]p -b Vl^/iS/cc - 10 “^+ ^JXtT+ci/\p, 

(C13a) 

p= CmP + w AtKcr-/4[e“*'^- c|),Cc-bH. C.,p] -br+X>[c|„]p 

-b (1 - Xt)T- Vlcmlp + 'D[^/K/,Cc - ie“^- i/AtT-c^jp, 

(C13b) 

respectively, with the definition ip± = aig(p±). If we 
additionally choose Aic = ±a;)[f the resonant terms have 
phases ip± = 0 and the resonant scattering rates are r± = 
4g^/Kic, while the off-resonant rates Tip = eT with e = 
l/[l-b(4wm/Kic)^]- Equations (C13) then lead to (5) and 
( 6 ). 


Appendix D: Evaluation of the correlation 
parameter S' 

To evaluate the quantity S in (1) for the bipartite sys¬ 
tem consisting of the two light pulses, we need to model 










the sequential measurement of the pulses using two mi¬ 
crowave cavities containing a qubit. The cavity modes 
then effectively constitute systems A and B. Assuming 
that initially the mechanical mode is in a thermal state 
with a mean occupation number uq and both cavities are 
in the vacuum state, this allows us to find the final state 
Pab of modes A and B before the qubit measurements. 
This is achieved by first integrating (5) for a duration ti 
and then (6) for a duration T 2 . The state pab is then 
used to evaluate S, while the mechanical mode is traced 


out. (In fact the mechanics should nearly factor out from 
the rest of the system at this point.) 

Under the dynamics described by the adiabatic mas¬ 
ter equations (5) and (6) the covariance matrix = 
-h XiXk) - {Xk){Xi) of the vector X = 
{xa,ya,Xb,yb) evolves as described by the differential 
Lyapunov equation S = FS -)- SF"'’ -|- N. For the blue- 
detuned pulse we find 


/ (1 6)rsq 

0 

0 

0 ^ 


f 

Tsq 

0 

y/ Vt^c^sq 

0 ) 

0 

(l sq 

0 

0 

, N= - 
’ 2 

0 

fsq 

0 

y/ 'Ht^c^sq 

sq 

0 

-Kc 

0 


' "x/ Vt^c^sq 

0 

Kc 

0 

V 0 

2\/ Vt^c^sq 

0 

-Kcj 


[ 

0 

\/ ‘Ht^c^sq 

0 

Kc j 


while for the red-detuned second pulse we have 

/-(I - e)rbs - 7 m 0 0 

p ^ 1 I 0 -(1 - e)rbs - 7m 0 

2 I 2 -^^ Pt^c^hs 0 

V 0 -2y/r]tKcTbs 0 

with fi = (1 -I- e)ri -I- 7m(2h -I- 1). To find S it suffices 
to find evaluate the covariance matrix at the end of the 
pulse sequence. 
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